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Abstract 

In this work we outline a program for lattice QCD that would provide a first step toward under- 
standing the strong and weak interactions of strange baryons. The study of hypernuclear physics 
has provided a significant amount of information regarding the structure and weak decays of light 
nuclei containing one or two A's, and S's. From a theoretical standpoint, little is known about the 
hyperon-nucleon interaction, which is required input for systematic calculations of hypernuclear 
structure. Furthermore, the long-standing discrepancies in the P— wave amplitudes for nonleptonic 
hyperon decays remain to be understood, and their resolution is central to a better understanding 
of the weak decays of hypernuclei. We present a framework that utilizes Liischer's finite-volume 
techniques in lattice QCD to extract the scattering length and effective range for AA^ scattering in 
both QCD and partially-quenched QCD. The effective theory describing the nonleptonic decays of 
hyperons using isospin symmetry alone, appropriate for lattice calculations, is constructed. 
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I. INTRODUCTION 



Nuclear physics is a fascinating field that remains poorly understood at a fundamental level. 
While the underlying Lagrange density describing all strong interactions is well established 
to be Quantum Chromo Dynamics (QCD), its solution in the nonperturbative regime leads 
to many unexpected and intriguing phenomena when the fundamental constants of nature 
take their physical values. Furthermore, one knows effectively nothing about how most 
strongly interacting systems behave as the fundamental constants are allowed to move away 
from their physical values. A puzzling feature of nuclear physics is the unnaturally large 
size of the S"— wave nucleon-nucleon scattering lengths; i.e. one would like to understand 
why low-energy QCD has chosen to be very near an infrared unstable fixed point of the 
renormalization group, as this requires a non-trivial conspiracy between the quark masses 
and the scale of QCD, Aqcd- While nuclei composed of neutrons and protons are well studied 
throughout most of the periodic table, and further efforts are proposed to study very short- 
lived nuclei that impact stellar environments, relatively little is understood about nuclei that 
contain one or more valence strange quarks: the hypernuclei. Understanding of the structure 
and decays of hypernuclei from first principles remains at a primitive level, and it is possible 
that the study of hypernuclei may help in unraveling fundamental issues regarding ordinary 
nuclei in terms of QCD. 

Since the A hyperon is the lightest among the hyperons (ttia ~ 1115 MeV), most theo- 
retical and experimental efforts have concentrated on bound systems of non-strange baryons 
and one A. These single-A hypernuclei can be experimentally produced by strangeness- 
exchange reactions N{K, 7c)A at CERN, BNL, KEK and Da$ne, by strangeness associated 
production reactions n(7r, K)A at BNL and KEK, or by the electroproduction mechanism 
p{e,e'K~^ )A at JLab (also planned at GSI for 2006) yj. These reactions usually populate 
highly-excited hypernuclear states, which decay electromagnetically to leave the hyperon in 
the lowest energy state, lsi/2. Once in the ground state, the system decays weakly, with- 
out conserving strangeness, parity or isospin. For very-light nuclei, such as ^H, the decay 
proceeds mainly through the mesonic decay mode, A — A^vr, which is analogous to the 
decay of the A hyperon in free space. These pionic decay channels are highly suppressed 
when the A is embedded in a nuclear medium, as the momentum of the outgoing nucleon 
(~ 100 MeV) is small compared to the typical Fermi momentum (~ 270 MeV for infinite 
nuclear matter). Non-mesonic (multi)nucleon induced decay channels, mainly AA^ — A^A^, 
are then responsible for the decay of the A hyperon in the medium, and become dominant 
for p-shell and heavier systems (note that even^ for ^He, the non-mesonic decay channel 
amounts to 40% — 50% of the total decay rate 2]). The lack of stable hyperon beams (the 
A lifetime is ta ~ 2.6 x 10~^° s) makes it difficult to perform a study of the strong and weak 
hyperon-nucleon (YN) interaction in free space. Therefore, hypernuclear structure and decay 
experiments presently provide the only quantitative information on the YN interaction. 

Experimental information on hypernuclear decay results from the measure of total life- 
times, partial decay rates (typically from the detection of a proton in the final state, and 
more recently from the detection of nn and np pairs in coincidence) and angular asym- 
metry in the emission of protons from the decay of polarized hypernuclei. However, the 
existing data have large error bars, in part due to poor statistics. On the theoretical side, 
one has to fold a particular model for the weak AA^ NN reaction into a many-body 
system. Realistic ground-state wave functions in the S = —1 sector have been derived 
for A < 4 hypernuclei using the Faddeev-Yakubovsky formalism 0], and up to j9— shell 
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hypernuclei using a three- and four-body cluster structure solved with the Gaussian-basis 
coupled-rearrangement-channel method The complete description of the final few-body 
scattering states in hypernuclear decay has only been performed for A = 3 using the Faddeev 
formalism j^. For heavier systems, finite nuclei are typically studied within a shell-model 
framework, while the final state is treated as a distorted wave function describing the two 
weakly-emitted nucleons moving away from a residual {A — 2) core. In any of those calcula- 
tions, a model for the YN interaction (typically a one-meson-exchange model) is assumed in 
order to describe the strong two-body dynamics. The development of Effective Field The- 
ories (EFT) in the strong S = sector, motivated recent model-independent work in the 
S = —1 sector. An EFT 0] was constructed at next-to-leading order (NLO) for momenta 
below pion-production threshold, and used to study scattering YN observables, as well as 
hyperon mass shifts in the medium, via a low-density expansion. The existence of the hy- 
pertriton was used to constrain the AA^ singlet- and triplet-scattering lengths. However, a 
description of bound states using this systematic two-body formalism does not yet exist. 

Regarding the \AS\ = 1 AA^ interaction, all the models developed to date agree in 
describing the long-range part of the weak transition through the exchange of Goldstone 
bosons, 71, K and t]^. For the short-range part of the interaction, some studies adopted 
the exchange of heavier mesons 0, 0] , the vector p, oj and K* mesons, while others used 
a quark-model type description [9|. Although some progress has been achieved during the 
last ten years, none of the existing models give a satisfactory explanation of the decay 
process. Construction of a model-independent EFT description of these decays was initiated 
in Refs. jiol . The study of the AS" = 1 AA^ interaction has the additional charm of 
allowing access to both the parity-conserving (PC) and the parity-violating (PV) part of 
the weak interaction, in contrast to its AS" = A^A^ partner, where the PC weak signal is 
masked by the strong PC amplitude which is 10 orders of magnitude larger. 

A major tool for strong-interaction physics that will play an ever-increasing role in nu- 
clear physics is lattice QCD. First-principles simulations on the lattice are the only known 
way to perform rigorous calculations of observables and processes where the QCD strong in- 
teraction contributes. Present-day calculations in lattice QCD are restricted to fairly large 
quark masses, relatively small volumes and relatively large lattice spacings. However, as 
time progresses the lattice-quark masses are decreasing, lattice sizes are increasing, and 
lattice spacings are decreasing. Furthermore, it has been appreciated that unphysical sim- 
ulations performed with partially-quenched QCD (PQQCD) allow for the determination of 
unquenched QCD observables, and can be performed sooner than fully unquenched simula- 
tions. While significant effort has been placed in calculating observables in the meson sector 
with impressive success, systems involving more than one baryon have largely been ignored 
with the exception of one pioneering effort to compute the nucleon-nucleon scattering length 
in quenched QCD 

In this work we establish the framework with which lattice (PQ)QCD can determine ba- 
sic strong-interaction parameters important for hypernuclei: the scattering parameters for 
AA^ AN. Given the presently large experimental uncertainties in this scattering ampli- 
tude, a lattice calculation will provide an important contribution to refining the potentials 
that are input into hypernuclear-structure calculations. Another component of this work 
is to explore the fundamental ingredients of calculations of weak decays of hyperons that 



The ry-exchange mechanism is omitted in some works due to the smahness (a factor of 10) of the rjNN 
coupHng constants. 
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must be understood before we are able to claim an understanding of the weak decays of 
hypernuclei themselves. Such an understanding has eluded theorists for decades; while one 
appears able to describe the S"— wave amplitudes in the hyperon nonleptonic weak decays, 
the P— wave amplitudes are very-poorly predicted. Lattice QCD is perhaps the only hope 
for gaining an understanding of these fundamentally-important weak processes. 



II. STRONG AN SCATTERING 

In predicting the structure of hypernuclei, the hyperon-nucleon potential or the hyperon- 
nucleon scattering amplitude is required input. A compilation of the world's experimental 
data on elastic hyperon-nucleon scattering can be found at the Nijmegen "NN-Online" web- 
site fl3]. The cross-sections in particular energy regions for these processes have large 
uncertainties associated with them as the total number of events is quite small ^. To get 
a sense of the uncertainty in the effective range parameters, consider that Ref. [14] quotes 
the following one-standard-deviation bounds for the 5*— wave scattering lengths and effective 
ranges: 

0.0 > a^'^") > -15 fm 0.0 > r^'^") > 15 fm 

-0.6 > a('^i) > -3.2 fm 2.5 > r^''^^) > 15 fm . (1) 

By comparison, a contemporary study jl^ finds (best fit) a^^^°'' = —1.8 fm, r^^^°^ = 2.8 fm, 
and a*^^'^!^ = —1.8 fm, r'^^'^^^ = 3.3 fm. Unfortunately, the meaning of these effective-range 
analyses is unclear. They arise from a four-parameter fit to data with very-poor statistics at 
energies where one would not necessarily expect shape parameters and other higher-order 
terms in the effective-range expansion to be small. Moreover, the effective-range parameters 
fit to the data are highly correlated. While the bounds of eq. are therefore suspect, they 
do encompass many model-dependent analyses of hyperfragments, production reactions and 
potential models [l6|. What then is known? Given the current state of experiment, it is 
probably safe to say: (i) there is no hyperdeuteron and therefore a^^^"'' < and a'^^'^^^ < 0. 
(ii) consistency of potential models with the spin of the hypertriton implies la^^*^"^! > la*-^"^^^!. 

With the lack of experimental effort in the direction of Lambda-nucleon scattering, lattice 
QCD may provide a rigorous, first-principles calculation of these cross-sections over some 
energy interval with smaller uncertainties than those determined experimentally in the near 
(and possibly far) future. In this section we put in place the theoretical framework required 
to connect possible future lattice calculations with data in the low-energy regime, below the 
pion-production threshold. 

The past decade has seen remarkable progress in the development of an EFT description 

of nucleon-nucleon (NN) scattering; ^Two quite different power-countings were developed 

for these processes, Weinberg H Q and KSW power-counting B H H Q, 

with neither providing a complete description of scattering in all partial waves. A third 



power-counting, BBSvK counting 25|, which is essentially a synthesis of these countings 
appears to be both formally consistent and provides a convergent expansion. The difficulty 
in establishing an EFT description is due to the large scattering lengths in both the ^5*0 
and ^5*1 channels, or equivalently, the near threshold (un)bound states in both channels. 



^ We are grateful to Rob Timmermans for discussions about the state of hyperon-nucleon scattering data. 
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Such states cannot be generated at any order in perturbation theory, and require the non- 
perturbative resummation of at least one operator in the theory. Physically speaking, such 
states result from a fine-tuning between kinetic and potential energy, and imply that the 
EFT is in the proximity of an unstable infrared fixed point jl^ (and thus implies that QCD 
has an infrared unstable fixed point). 

There is no reason to suppose that the hyperon-nucleon sector is near an infrared fixed 
point. That is to say, while the bounds of eq. ((H) do not preclude unnaturally large scat- 
tering lengths, we would expect the scattering lengths to be of natural size, and for there 
to be no AA^ bound states near threshold. Therefore, one expects the EFT that describes 
the interaction between hyperons and nucleons to have a power-counting based on the en- 
gineering dimensions of the operators that appear in the effective Lagrange density, and 
hence to be precisely the power-counting of Weinberg 17, ^- As a well-defined procedure 
exists, due to Liischer 27, li^, for extracting the low-energy elastic scattering parameters 
from finite-volume lattice QCD calculations, we set-up the low-energy EFT describing AA^ 
scattering in Weinberg power-counting^. Given the Liischer framework, we use the EFT to 
compute the scattering length and effective range for elastic AA^ scattering to NLO in the 
EFT expansion ^. These parameters will be sufficient to describe AA^ scattering for energies 
below that necessary for inelastic processes AA^ — > SA^ and AA^ — * AA^vr. (See also Ref. p.) 
Our results contain the leading light-quark mass dependence of low-energy AA^ scattering 
and can be used to extrapolate from lattice-quark mass values to nature. 



A. AA^ ^ AA^ in QCD 

If one could perform a lattice QCD simulation at the physical values of the quark masses, 
then, barring the advent of new experiments, there would be no need to construct the EFT 
for the AA^ AN S-matrix at low energies. However, near-future lattice simulations will 
be partially-quenched with unphysically large values of the sea-quark masses. Therefore, at 
least in the short term, the EFT construction will be needed to extrapolate to QCD. 

In defining the nonrelativistic fields appropriate for computing the strong interactions 
between the A and the nucleon we remove the nucleon classical trajectory from the nucleon 
field, and the A classical trajectory from the A field ^. We can perform these redefinitions 
as we will work with an SU{2)l ^ SU{2)b. chirally-invariant theory, as opposed to its three- 
flavor analog. The Lagrange density describing the free-field dynamics of the nucleon, A and 
S is 



^ Weinberg's power-counting does not, in general, require resummation of interactions. 

^ Note, however, that we power-count the amphtude directly, rather than the potential. 

^ The three-body system with one A-hyperon has been investigated in an EFT without pions |22j. 

^ Usually, one removes only the classical trajectory associated with one of the heavy-hadron fields. However, 

we will work with two-flavor xPT in which strange-baryon number is conserved at each interaction, and 
hence the AA'^ residual mass term, Aaat can be removed from the theory by a phase redefinition, leaving 
only the SA residual mass term. 
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where A is an iso-singlet while the nucleoli iso-doublet and the S iso-triplet are defined as 

J]V^) ■ (3) 

As we are interested only in AA^ scattering up to NLO, we need not consider the Weinberg- 
Tomazawa type interactions. The complete kinetic energy terms for the A^ and S have 
additional multi-pion interactions, but as the A is an isosinglet its chirally-invariant kinetic 
energy term is given in eq. The leading-order (LO) interactions between the baryons 
and the pions are given by (to linear order in the pion field) 

C = ^ N^a- VMN + ^ ( A"^Tr [ a ■ VME ] + h.c 

+ ^ Tr [ St [a ■ VM , S ]] , (4) 

where / ~ 132 MeV is the pion decay constant, qa ~ 1-27 is the nucleon axial coupling 
constant, (^as ~ 0.60 is the AS axial coupling constant, while g^.^. is the SS axial coupling 
that cannot be determined directly from experiment and which we will discuss later. M is 
the matrix of pions given by 

M - f'l^ _X,^) . (5) 



TT — TT 



At LO in the expansion, there will be contributions to 5— wave scattering from four- 
baryon operators of the form 

C = AACf ( A^ P('^») A^)^ ( A^ P('^°) n) + AACf ^) ( A^ P('^^) Nf ( A^ p('^^) n){6) 

where P^^^o) = ^ ^^-^ an^l p('^'^i) = ia2(y"' are spin-projectors for the ^5*0 and ^5*1 channels 

respectively, and the coefficients aaCq ' be determined. At next-to-leading order 

(NLO) there are also contributions from four-baryon operators that involve a S baryon, 

C = sACf ( S^(^r2) P(^^«) Nf ( A^itn) P(^^«) n) 

+ EAC^r^ ( ^^{^T2) P^''^^ Ny ( A^izr,) P'-''^^ n) + h.c. , (7) 

where the coefficients t,aCq^°' are to be determined. 

We have recognized and implemented a particular hierarchy of mass-scales in order to 
simplify our expressions for the scattering length and effective range. By returning briefly 
to SU{3), we recall that the SA mass sphtting, Aas ~ rrig is proportional to the light- 
quark masses, in particular, the strange quark mass, and numerically is Aas ~ 74 MeV. 
In contrast, the pion mass is ~ and numerically is m,r ~ 140 MeV. Therefore, in 
what follows we will neglect Aae compared to m^. However, we will retain terms of the 
form A/2yUAAfAAs (where fi\N is the reduced mass of the AA^ system) compared to m^, as 
they are of the same order in the SU{3) chiral expansion. While we can (somewhat) justify 
this expansion from a formal point of view, one should keep in mind that the expansion is in 
factors of ~ 2. Of course the sole purpose of this expansion is to provide simple closed- form 
expressions for the effective-range parameters; the exact numerical expressions are easily 
obtained and may in fact be necessary in certain regions of the Co parameter space. 
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FIG. 1: Diagrams that contribute to AN scattering at LO (top) and NLO in the EFT. The shaded 
blobs are vertices from eq. 



A straightforward calculation of the diagrams in Fig. ^yields the scattering length a^^^°'^ 
and effective range "^"^ in the ^5*0 channel: 



aaC'o 



{'So) 



4^ V^^^° 



('So) 



^{i5o) S^-sA gA /iATV v"^ + V^T. + ml 

t^n : — TT^ 



1] + mT, 



2(?7 + m^)2 



1 



277 



aaC'o 



57r 



(^SAL-o 



/iAjV 



saC'o 



27r/2 



+ 



12(r/ + m^)4 



where ?7 = v^2yUAArAAs. We see that there are two unknown constants, aaCq ° and saC'o "^"^ 
and two observables, a^^^^"^ and r^'^^''^ In the spin-triplet channel we find 



6(77 + m^) 



.('So) 



(8) 



aaC'o 



47r 



2nP 



r] + 



^a'ik 9a f^AN 2?7^ + Ar]^m^ + 6r]ml + 3m^ 



47r/4 



2{r] + m^)2 



(9) 



and 



/iAAfTT 



27r 



AA<-^0 



2 r 



3 / 
8^ 



fl'SA gA fJ'AN 377^ + Or^m^ + 8ml 

^On r 77^ TT^ 



?7 ° 27r/2 6{ri + m^)^ 



3fl'lA Z^AAf 46?]^ + gir^^m^ + AArjml + llm^ 



47r/4 
see t 

observables, a'-^'^^-' and r*^^'^^^. 



36(r/ + m^)4 

where we again see that there are two unknown constants, aaC'o and saC'o , and two 



(10) 
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It is worth commenting on the viabihty of the power-counting scheme that we have 
chosen for the AN system. Consider the ^Sq effective-range parameters of eq. (jH)). If we 

vary the Cq^°'' coefficients over natural values (i.e. ~ 1/ p), then there are large regions of 
parameter space where there is poor convergence in the scattering length, and the effective 
range is anomalously small and negative. This may indicate that a next-to-NLO (NNLO) 
calculation is required, which includes, inter alia, contact operators with two derivatives 
(the C2 operators) and contact operators with an insertion of the quark-mass matrix (the 
D2 operators). However, it may also be the case that we should apply our power-counting 
scheme to a potential, rather than the amplitude directly, and resum using the Schrodinger 
equation, as in Weinberg power-counting in the NN system. The appropriate power-counting 
will ultimately have to be decided in conjunction with lattice QCD simulations of the AA^ 
system. 



B. KN AN in PQQCD 

While one always hopes for fuUy-unquenched lattice simulations at the physical values of 
the light-quark masses, it will always be the case that partially-quenched simulations of 
the observables of interest will be performed first. This is for the simple reason that the 
partially-quenched simulations, in which the masses of the sea-quarks are larger than those 
of the valence quarks, take less time than the fuUy-unquenched simulations. During the 
past year or so a significant number of single-nucleon observables have been explored in 



partially-quenched chiral perturbation theory (PQxPT) j30|, |3l|, |3^ with an eye to making 
connections between partially-quenched lattice simulations and nature. The EFT describing 
NN interactions has also been partially quenched ^33. .34] . 

The partial quenching of the AN scatte ring amplitude has features from both the 
partially-quenched NN scattering amplitude |3 and from the AqAq potential [3^. The 
local four-baryon operators in PQQCD are the same as those in QCD because ghost-quark 
and sea-quark number are conserved. It is for this reason that the pion box-diagram is the 
same in PQQCD and QCD. The Crossed-Box diagram is modified by partial quenching, but 
this effect enters at NNLO, one order beyond the order to which we are working. Unlike the 
AqAq potential |3^, there is a contribution from one-hairpin exchange (OHPE) at LO and 
also from one- and two-hairpin exchanges at NLO, as shown in Fig. |21 

* N 



N> 



^ A A 
A A f 



FIG. 2: The LO and NLO contributions to AN scattering from hairpin-exchange in PQxPT. The 
shaded blobs are vertices from eq. 111^) . 



The SU{2)l SU{2)r chiral symmetry of two-flavor QCD is extended to the graded 
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SU(4\2)l ® SU{4\2)ii chiral symmetry of PQQCD, as discussed extensively in Refs. jsi], 
Is^ . In this theory, in addition to the hght pions, hght ghost-mesons and sea- mesons that 
contribute in diagrams with single-pole propagators, there are also contributions from rj- 
exchange with a double-pole propagator of the form 



mz 



(g2 _ 777,2 _^ ^(r)2 



(11) 



where m^^ is the mass of a meson composed only of sea quarks. We are working in the 
isospin limit where the sea-quarks are degenerate and the valence quarks are degenerate. 
This double-pole arises because the singlet of SU{2) is not the singlet of 5'f/(4|2) [39]. At 
LO in the chiral expansion, the interaction between the rj and the baryons is described by a 
Lagrange density of the form 



(12) 



where the axial couplings, Qq^'^'' are to be determined from lattice simulations. The analo- 
gous contributions to the AqAq potential js^ are absent because the axial coupling of the 
hairpin to the Ag vanishes due to heavy-quark symmetry, and is non-zero only at order 

l/mq [i^. However, we expect the coupling g^^^^ to be significantly smaller than qy.k based 
on large- A^c arguments. In the large- A^p limit, the "size" of the matrix element of an op- 
erator with isospin / and spin J is of order 1/N]i~'^^ 37, ll^. Hence we expect this axial 
coupling, which has / = and J = 1, to be suppressed by a factor of ~ 3 compared to (7sa, 
which gives it a natural size of (/q^^ ~ 0.2. 

We find that the leading contributions to the scattering length and effective range in the 

channel due to partial quenching are 



' 271 
1 



(A) (N) 



ml, - ml 



1287r/^ 



271 



AA'-^O 



/iAATTT 

and in the ^Si channel are 



2 (A) (N) 2 



— mz 



1 

yUAAfTT 



mz 



(13) 



— aaC'o 



27T 
LAA<^0 



) 96 % f^AN m^^ - m^ 



+ 



m^ 



1287r/^ 



mi 



2 (A) (N) 2 2 

9o 9o m^s - 



3f 



ml 



(14) 



It is important to note that the additional contributions to the scattering lengths arise from 
the NLO diagrams in Fig. |2l and as required, there is no contribution from OHPE (due to 
the derivative coupling). By contrast, the contribution to the effective ranges begins at LO 
from OHPE. Therefore, partial-quenching will perturbatively modify the scattering lengths, 
but overwhelmingly dominate the effective ranges, except very near the QCD limit. Given 
that extractions of both the scattering length and effective range are required to determine 
the four-baryon, momentum-independent leading-order interactions, partial-quenching may 
provide a serious obstacle to this procedure. 
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C. AA^ AN at Finite Volume 

The Maiani- Testa theorem ji^ precludes determination of scattering amphtudes from 
Euchdean-space Green functions at infinite volume away from kinematic thresholds. While 
appearing to be a significant setback for lattice QCD for which only Euclidean-space Green 
functions can be computed, it was realized by Liischer ji^, 113 ^^^^ '^^^ access 2 —>■ 2 
scattering amplitudes from lattice calculations performed at finite volume. Liischer's work 
generalizes a result from nonrelativistic quantum mechanics to quantum field theory. 
The results are obtained by assuming that the lattice spacing, 6 '', is very small and the 
lattice size, L, is much larger than the range of the potential between the two particles. 

During the past few years, the pionless effective field theory, EFT(/) 11113, description 
of NN scattering at momenta much less than the pion mass has been studied extensively, 
and it can be applied to describe the finite- volume construction described above ^. In 
EFT(^) the interactions between nucleons are described by a series of delta-functions and 
their derivatives. While the two-nucleon sector reproduces identically effective-range theory, 
when gauge fields are included one has a systematic and calculable scheme with which to 
determine all low-energy multi-nucleon processes. This theory has been successfully applied 
to quite a wide variety of two-nucleon and three-nucleon processes^. In order to establish a 
power-counting at the level of diagrams in systems with unnaturally-large scattering lengths, 
the power-divergence subtraction scheme (PDS) was introduced to deal with the power-law 
ultra-violet divergences that appear in the non-relativistic theory at the one-loop level. It is 
somewhat useful to see how one recovers Liischer's relations between the energy-levels of the 
quantum system in a periodic box and the scattering parameters of the continuum theory in 
EFT(^), and this is presented in Appendix I. As re-derived in Appendix I, the energy of the 
lowest-lying scattering state of the two-particle system with zero center of mass momentum 
in a periodic box of sides L that is much greater than any of the scattering parameters, such 
as the scattering length and effective range, is ji^, lisj ] 



where the coefficients ci,2 are ci = —2.837297 and C2 = 6.375183, and where we have used 
the conventional definition of scattering length and effective range 



The energy of the next level which transforms in the Ai representation of the cubic group 



^ In this work we use b instead of a to denote the lattice spacing. Clearly this will be considered a sacrilege 
in the lattice community, but denoting the lattice spacing by the traditional a would lead to confusion 
with the scattering lengths which are also traditionally denoted by a. 

^ This discussion will be the same as that for pseudo-potentials |4J| as the "physics" is the same. However 
the EFT(/) construction allows for generalization to processes including gauge fields. 

^ One may worry about potential violations of unitarity in EFT(/) for PQQCD. However, as pointed out 
in Ref. [S^ i one can use the methods of Ref . [s^ to show that all effects of partial-quenching are in the 
coefficients of the contact operators. Therefore, perhaps somewhat surprisingly, EFT('^) is unitary in 




(15) 



pcotS = 1 — rp^ -|- ... 

a 2 



(16) 



PQQCD. 



10 



is (U 111 



E 



27r^ 6tan5o 



( 1 + c;tan5o + C2tan2 5o + ...)+ 0{L-^) , (17) 



1 



where c[ = —0.061367 and C2 = —0.354156, and where 6q is the 5— wave phase shift evaluated 
at the unperturbed lattice momentum |p| = 2tt/L. It is important to stress that the effective- 
range parameters may be comparable to or greater than realistic box sizes at some values of 
the unphysical lattice quark masses. Eqs. (fT5|) and (fTTjl are then no longer applicable and 
one must use the exact formula for the energy levels, eq. (j38p . given in Appendix I. 

III. NONLEPTONIC WEAK DECAYS AND WEAK HYPERNUCLEAR DECAYS 

The weak decays of hypernuclei result not only from the one-body nonleptonic decays 
(mesonic) of the hyperon(s), e.g. A — > A^vr, that exists in free-space but also from two- 
and higher-body decays (non-mesonic), such as AA^ NN and ANN NNN. Naive 
expectations suggest that the one-body decays should dominate, however the non-mesonic 
processes are found to be comparable to the mesonic decays in nuclei, in part due to Pauli- 
blocking of the final-state nucleon in mesonic decays. 

The nonleptonic weak decays of hyperons in free-space have been very-well studied dur- 
ing the past many decades. While the 5*- and P— wave amplitudes that contribute to each 
process are very-well known experimentally, the P— wave amplitudes are still not well de- 
scribed theoretically. In order for one to claim that the weak decays of hypernuclei are 
theoretically understood we must first understand the nonleptonic weak decays of hyperons 
in free-space, as this underlies one of the two competing decay processes. Furthermore, the 
EFT that describes these processes is more than likely entangled in the theory that describes 
the higher-body amplitudes, and therefore is central to the description of all decay processes. 

A. Nonleptonic Weak Decays of Hyperons in Free-Space 

As the nonleptonic decays A —>■ pn^ , rnr^ , S+ nn^, —>■ Avr" and their isospin 

partners have been well studied for an extensive period of time, the experimental data and 
theoretical efforts to understand the data can be found in many textbooks, e.g. Ref. jisj ]. 
Therefore, we do not attempt an exhaustive discussion of the existing literature in this 
work. Angular-momentum considerations are sufficient to write the matrix element for the 
nonleptonic decay B — > B'n, Mbb'-k, using the notation of Jenkins ji^, as 



where k is the four-vector of the outgoing pion normalized to the pion three-momentum 
in the baryon rest frame. The factor of that appears in the coefficient in eq. ()18|) is 
the physical mass of the charged pion, and does not indicate that the amplitude vanishes 
in the chiral limit. The amplitudes A^^^ and A^^^ are the S— and P— wave amplitudes, 
respectively. They are expected to be the same order in the chiral expansion, and indeed 
experimentally, their magnitudes are of order unity. 



M = GFml+ Ub' \A^^^ + 2 S^'k^ A^^^ 1 U, 



B 



(18) 
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1. Nonleptonic Weak Decays of Hyperons in SU{3) 



Conventionally, the theoretical analysis of these decays is performed about the limit of SU (3) 
flavor symmetry, assuming the AI = | rule, although a more general analysis shows that 
the AI = I rule is well respected in these decays j4^. At LO in the chiral expansion, the 
weak decays are described by the Lagrange density 



L = GFml+ f [ hn Tr [B{h^, B} 



+ hp Tr 



B[h,B] 



(19) 



where 



h ^ with h 




(20) 



and where h transforms as a (8, 1), /i — ^ LhL^ , under SU{3)l ® SU{3)ji chiral transforma- 
tions. Expanding this interaction to linear order in the pion field, it is straightforward to 
show that the 5— amplitudes arising from the diagram in Fig. IHlor its analog are |4fil. l47l| 




FIG. 3: Tree-level diagram that gives the leading-order contribution to the S—wave amplitudes for 
the nonleptonic decay A — > pir^ . The solid square denotes an insertion of the leading-order weak 
interaction from eq. Iliy\) . 



A. 



-Att" 



hn — hp 

75 (3/iF - ho) 



(21) 



and the P— wave amplitudes arising from the diagrams in Fig. E] and their analogs are 



A Z 



A 



n P 



\ 71 
\ 



\ 

\ 

71 X 



FIG. 4: Tree-level diagrams that give the leading-order contribution to the P—wave amplitudes for 
the nonleptonic decay A pir^ . The solid-square denotes an insertion of the leading-order weak 
interaction from eq. I\iy\) while the shaded circles denote axial couplings. 
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(P) 



DihD-hp) (^D + F){hD + 3hp) 



Mn - Mso 3(M„ - Ma) 
D{hD + ^hF) F{hF-hD) {D + F){hD-hi 



3(M„ - Ma) M„ - Meo Me+ - Mj 



(5) ^ (D-F)(3/ip-M bD{hn + hF)_ .221 
y6(MA - Mho) V 3 Mh- - Ms- ' ^ ^ 

where D and F describe the matrix elements of the octet axial current operators between 
octet baryons. 

In table |l] we present the numerical values for the LO tree-level amplitudes, along with 
their experimental values. We choose to determine ho and hp from the 5*— wave amplitudes, 
which give central values hp, = 0.58 and hp = —1.40, and to use the tree-level values for the 
axial matrix elements, D = 0.80 and F = 0.50 Also, in table|l]we present the amplitudes 



TABLE I: Weak Amplitudes in SU{3) xPT at LO and NLO 



Decay 


^(^) LO 


^(^) NLO [50] 


^(^) Expt 


^(^) LO 


^(^) NLO [50] 


^(^) Expt 


A pTT~ 


1.48 


1.44 


1.42 ±0.01 


0.59 


-0.73 ±0.18 


0.52 ±0.02 


mT~ 


1.98 


1.89 


1.88 ±0.01 


-0.30 


0.46 ±0.21 


-0.06 ±0.01 




0.0 


0.01 


0.06 ±0.01 


0.16 


-0.18 ±0.21 


1.81 ±0.01 


Avr" 


-1.95 


-2.01 


-1.98 ±0.01 


-0.19 


0.52 ±0.29 


0.48 ± 0.02 



that are found at one-loop order (NLO) in the three-flavor chiral expansion (see also 
Ref. 0). 

The numerical values in table U make it clear that the theoretical analysis does not 
come close to predicting the P— wave amplitudes. Moreover, the one-loop computation in 



three-flavor xPT |47|, |49|, |5 



5ll | does not improve the situation and, arguably, makes the 
situation worse. There are a variety of reasons one can conjecture as a possible explanation 
of these discrepancies, e.g. SU{3) breaking is abnormally large, there are strong final state 
interactions, there are narrow hadronic resonances (e.g. pentaquarks) in some channels, and 
so forth. From a formal standpoint, it is possible that the three- flavor xPT expansion simply 
does not behave well for these observables, given that it is the kaon mass that is governing 
the convergence. Jenkins |4^ has argued that convergence is poor, not due to failure of the 
chiral expansion, but due to large cancellations at LO. In order to explore this possibility 
further, a NNLO calculation is ultimately required. As parameters that enter into the NLO 
calculations were determined from a global fit it is somewhat difficult to determine if the 
problems in the P— waves lie in only one of the decay modes or if the theory fails equally 
badly in all modes. 



1° Ref. |43 finds D = 0.79 ± 0.10 and F 
the available semileptonic decay data. 



0.47 ± 0.07 from a global fit of tree-level axial matrix elements to 
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2. Nonleptonic Weak Decays of Hyperons in SU{2) 



Given the absence of complete theoretical control in the nonleptonic decays in three-flavor 
xPT, we de-scope our efforts to two-flavor xPT with the hope of gaining some insight into 
which processes can be identified as problematic and which are under control, without having 
to consider the issue of SU{3) breaking. We restrict ourselves to the decays of the A and S 
as we are primarily interested in the study of hypernuclei containing these hyperons. 
At LO in the two-flavor chiral expansion, the weak Lagrange density is 

C = Gpml+f [ hA (m^) A + hs N^h^ + h.c. ] , (23) 

where 

= S}h and h = , (24) 

and where h, which transforms as h ^ Lh under SU{2)l SU{2)ii chiral transformations, 
is a (2, 1). As in the three-flavor case, we have assumed a AI = | weak interaction only. 

The weak interactions described in eq. ()23|) along with the strong interactions between 
the pions and the baryons in eq. (jH) give S"— wave amplitudes 

<l- = hA , A^'l^^. = , A^'l^, = , (25) 
and P— wave amplitudes 

.(P) _ 9AhA gA^hj^ .(p) _ 



"-P^- Ma - M„ Mp - Ms+ ' M„ - Ma M„ - M^o ' 

.{P) _ 9AT,hA gT,T,hj: QAhj: , . 

s+n^+ - m^-Ma Mn- Mso Me+ - Mp ' ^ ' 

We choose to fit the weak couplings /ia.e to the central values of the 5— wave amplitudes, 
and find that Ha = 1-42 and /is = 1.88. 

While the axial matrix elements qa = 1-27 and qat. = 0.60 are well-known experimentally, 
only an upper limit currently exists for S~ — >■ TPeV, and thus there is no direct measurement 
of gY.T,- Given the success of 5'[/(3) symmetry in relating the axial-current matrix elements, 
it is worth exploring its predictions for (yf^s, along with the best estimates for leading- 
order SU{?>) breaking. In the limit of exact SU{?>), the axial- coupling would be (/es = 
F, and thus = 0.5 at leading order. We focus on two quite different approaches to 
estimating the SU{?>) breaking. In one method the one-loop corrections to all the axial 
matrix elements were computed in three-flavor xPT including the decuplet as dynamical 
flelds, and parameters in the one-loop contributions were varied over a reasonable range 
in order to estimate the size of S't/(3) breaking effects. This analysis yields a range for 
5'ss, 0.35< (yfES^ 0.55. However, the analysis is somewhat ad hoc and depends upon the 
convergence of the three-flavor chiral expansion, but formally captures the leading-order 
SU{?>) breaking in the chiral limit. The other method js^, lE^I we use to estimate (?ss is 
to perturb about the \aige—Nc limit of QCD, a program that has been very successful in 
understanding known baryon properties, and has provided insight into the successes of the 
naive constituent quark model. This large— A'^c analysis of the axial matrix elements leads 
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to a range 0.30< (7ss^ 0.36 Given that these distinct estimates of (^ss have only a 
small overlap and neither are obviously incorrect, we take 0.30< Qt^t.^ 0.55 as the range of 
this axial coupling constant, and it is clear that it is quite uncertain at present. 

Inserting numerical values for the couplings and masses into the S*— wave and P— wave 
amplitudes in eq. and eq. (pUjl gives the numerical values for the nonleptonic amplitudes 
shown in table |nj At leading order, the P— wave amplitude in A decays is well predicted. 



TABLE II: Weak Amplitudes in SU{2) xPT at LO (c/ss = 0.30 ^ 0.55) 



Decay 


Theory 


^(•5) Expt 


Theory 


Expt 


A pn^ 


1.42 (input) 


1.42 ±0.01 


0.56 


0.52 ±0.02 


Ti~ — > riTT^ 


1.88 (input) 


1.88 ±0.01 


-0.50 ^ -0.14 


-0.06 ±0.01 


Ti~^ — > nvr"*" 


0.0 


0.06 ±0.01 


±0.42 ^ ±0.08 


1.81 ±0.01 



as it is in the three-flavor theory. However, we do not expect significant modifications to 
this result from higher orders in the two-flavor theory. Further, we see that the P— wave 
amplitude for S" nn^ is quite sensitive to the value of g^-^, which is presently quite 
uncertain. At the upper limit of the allowed range for g-^^:, the P— wave amplitude is close 
to what is observed. Finally, we see that the P— wave amplitude for S"*" — >■ nvr^ is not close 
to the experimental value for any reasonable value of gY,s, and theory underestimates the 
experimental value by ~ 4 in the best case. 



By extending the contents of table II in Ref. [s^ the coupHng gss is found to be 

2 



2 . 4 



(27) 



Usin g th e central values of the global fit values for the parameters a, b, C3 and C4 found in table Vll of 
Ref. |53| we arrive at the range for .gss- We thank Aneesh Manohar and Elizabeth Jenkins for assisting 
with this determination. 

For completeness we have also determined the axial coupling g-£ that appears in the interaction Lagrange 
density 



(28) 



where s denotes the iso-doublet containing the ^'^ and ^ . By extending the contents of table II in 
Ref. [5^ the coupling gs is found to be 



9= 



1,4 8 

—a — H — cq C4 

3 3^3 



(29) 



Usin g th e central values of the global-fit values for the parameters a, b^c^ and C4 found in table VII of 
Ref. [53 we find 0.26^ 0.30. Further, the one-loop analysis of Ref. |48j gives 0.18 5- q p,^, 0.36. These 
two estimates are compatible with each other, with, as before, the estimate from Ref. |48| having a larger 
range than that from Ref. 



15 



B. Nonleptonic Weak Decays of Hyperons and Lattice QCD 



Given the failure of both SU (2) and SU (3) xPT to describe the P— wave amplitudes of the 
nonleptonic decays, one would like to understand where the problems lie. Until this issue is 
resolved, one is unable to claim theoretical control in the nonleptonic decays of hypernuclei. 
Lattice QCD is, of course, the only rigorous theoretical technique that exists to compute 
these decays, and it is important to use lattice simulations as a diagnostic of the problem(s) 
with the xPT calculations. 

One should first compute the two-point functions that arise from the leading-order in- 
teractions in eq. ()23|) . The lattice calculation of this process injects energy at the weak 
four-quark operator in order to allow both the initial and final states to be on mass-shell. 
Hence one must consider additional operators in the EFT that are usually discarded as sur- 
face terms. In the SU{2) expansion, the mass difference between the A or E and the nucleon 
is taken to be a scale that is fixed in the chiral limit, with dependence on the quark masses 
arising at higher orders in the chiral expansion, i.e. the strange quark mass is fixed. One 
therefore has the free Lagrange density 



£ = Nkv-DN + AkvdA - Aa^tAU + 



(30) 



At NLO there are contributions from 



Gprnl^ [ ) IV d[ (Nh^) A ] + jf Aa^ (Nh^) A 



, -(i) ■ a 
+ jj^ IV ■ a 



+ ^ AsTvA^S/i^ + h.c. 



(31) 



where the surface terms vanish when energy and momentum are conserved at the effective 
weak vertex. However, in lattice calculations, where both the hyperon and the nucleon will 
be on mass-shell, these operators make a non-zero contribution to the amplitude at tree- 
level. This contribution must be included when attempting to extract the weak couplings 
^A,s from lattice calculations. Therefore, including such terms, the weak matrix element 
between a A and a nucleon arising from the diagrams shown in Fig. El becomes 



A 



FIG. 5: Tree-level diagrams that contribute to the weak A ^ N matrix element, including the 
leading-order surface term. The dark shaded square denotes an insertion of the weak vertex, h\ 
from eq. \2'^) . while the light shaded square denotes an insertion of the weak vertices, ' ; from 
eq. h'Jl)) . The crossed- square denotes an injection of energy at the weak vertex. 



{N (p 



,■(1) 



/:weak |A (p = 0)) = Gprnl^ [f hA + A^N (j? 



+ if 



UnUa- (32) 



Given that the P— wave amplitude for A — > pvr" is in good agreement with data, one 
expects that the value of extracted from the lattice will be in agreement with that de- 
termined from the 5*— wave amplitudes. However, it is not clear what one will find for /i^ 
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given the significant uncertainty in g-^^ and the complete failure to reproduce the P— wave 
amplitude for S+ — > mr^. Therefore, a lattice determination of /is is highly desirable. 
Furthermore, given that there is a significant uncertainty in ^f^s and this uncertainty propa- 
gates through into a large uncertainty in the P— wave amplitudes for the S decays, a lattice 
determination of g^-E is also highly desirable. 

A direct lattice measurement of both the S"— wave and P— wave amplitudes would also 
be very important. Assuming that the weak couplings /ia,s and the axial couplings are all 
consistent with the above discussion, then we need to understand the origin of the failure 
of the P— wave amplitudes. Given the progress in extracting weak amplitudes for K tttt 
using the finite-volume techniques described by Liischer and Lellouch [5^ one can hope for 
a direct lattice determination of these amplitudes at some point in the future. While the 
actual value of the amplitude should agree with experiment one should be able to disentangle 
the various contributions to the amplitudes and understand where the failure originates in 
the low-energy EFT. 



C. Weak Decays of Hypernuclei 

A systematic study of the decays of hypernuclei provides an even greater challenge than 
describing the nonleptonic decays of hyperons in free-space. It will certainly be the case 
that the interactions that are ultimately needed to describe the nonleptonic decays will also 
provide a contribution, and likely a leading-order contribution, to the multi-baryon weak 
interactions such as AA^ NN. For instance, the diagram shown in Fig. IHl arises from 




FIG. 6: An example of a contribution to An — > nn. 



inserting the weak two-point amplitude on an external leg of the strong nn nn scattering 
amplitude. Unfortunately, there is presently no consistent EFT that is applicable in the 
kinematic regime for AA^ NN where the energy release for an initial AA^ state at rest is 
E ~ 176 MeV. If an EFT exists for the strong processes, such as AA^ AN and A^A^ NN 
at these relatively high energies, one hopes that the higher-dimension weak operators are 
subleading compared to the two-point function, but given the uncertainty in the nonleptonic 
hyperon decays and the absence of a strong-interaction EFT, it appears premature to make 
any such power-counting arguments at this stage. 

However, having made all of these negative statements based on our current theoretical 
understanding, there has been an important phenomenological analysis 0, 11 1 recently 



completed by Parreno, Bennhold and Holstein (PBH) that gives one hope that an EFT can 
in fact be constructed to describe these processes. In PBH, contributions to the weak decay 
from weak-one-pion exchange (WOPE) and weak-one-kaon exchange (WOKE) were included 
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along with contributions from local As = 1 weak ANNN interactions. The power-counting 
used in determining the order of the various contributions is based on Weinberg's power- 
counting of NN interactions in which the number of powers of the small momenta is directly 
related to the order in the expansion. The weak interaction computed to a given order 
in Weinberg power-counting is then dressed with NN strong interactions in the final state 
and the initial-state interactions are taken into account by a correlated AA^ wavefunction 
determined from Nijmegen models. 

The usual divergence problems with local four-baryon interactions arise in this calculation 
and are dealt with in PBH by extending the interaction in coordinate space to a Gaussian 
form with the same spatial integral (of course, the dependence upon the form of the regulator 
must be, and is, higher order in the power-counting). The results of a global fit, and data 
can be found in table IIIIl 



TABLE III: Weak decay observables for ^He, j^^B and taken from Ref. l3, U |- The theoretical 
numbers correspond to including WOPE, WOKE along with the leading-order parity- violating and 
parity-conserving ANNN interactions with the short-distance parameters globally fit to the data. 



Observable 


Theory (best global fit) [10, U} 


Expt 


r(iHe) 


0.44 


0.41 ±0.14 , 0.50 ±0.07 


n/p (iHe) 


0.55 


0.93 ±0.55 , 0.50 ±0.10 


^(iHe) 


0.24 


0.24 ±0.22 




0.88 


0.95 ±0.14 


n/p (i^B) 


0.92 


-■-■^^-0.48 


^a^B) 


0.09 


-0.20 ±0.10 


ra^c) 


0.93 


1.14 ± 0.2 , 0.89 ± 0.15 , 0.83 ± 0.11 


n/p (i^C) 


0.77 


0.87 ±0.23 




0.03 


-0.01 ±0.10 



It is found that the overall goodness of the fit encapsulated in the is independent 
of the range of the short- distance Gaussian regulator (over the range explored) and also 
independent of the strong-interaction potential used to compute the final-state interactions. 
This is very encouraging. The numerical values of the bare coefficients in the weak potential 
depend upon the range of the Gaussian regulator, as one would expect, but the dependence 
of the observables is minimal. Therefore, it appears that a systematic EFT description of 
the weak decays should be possible in the near future, and it may be as simple as Weinberg 
power-counting in many of the decay channels. Clearly, this exciting possibility requires 
further exploration. 

It goes without saying that it would be very nice to see a calculation of the weak AN — >■ 
A^A^ amplitude directly from lattice QCD (a 5-point function!). While one knows how to 
extract the weak decay amplitudes for K ^ titi from the lattice, part of which is tuning 
the lattice volume so that one of the vrvr energy states on the lattice is degenerate with a 
kaon at rest, no such formulation currently exists for extracting weak amplitudes for 2 — >■ 2 
processes. In principle we do not see any formal road blocks to constructing a framework 
analogous to that of Liischer and Lellouch, however, given that the one-body decays are not 
yet understood we do not see any urgency in developing such a framework. 
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IV. CONCLUSIONS 



We have explored some aspects of strong-interaction physics that impact our understanding 
of the structure and decays of hypernuclei and that could be addressed by lattice QCD 
simulations in the near future. In the few-nucleon sector, lattice QCD will provide an 
understanding of the strong dynamics, such as the fine-tuning that exists in the two-nucleon 
system, but it will not improve the precision of input into nuclear calculations of many- 
body observables. However, in the hypernuclear sector where experimental measurements 
are scarce and quite uncertain, lattice QCD may be in a position to compete with and 
ultimately surpass experiment and thus greatly aid progress in the field. Wc have discussed 
quantities in both the strong and weak sectors whose calculation in lattice QCD would have 
great significance. 

To summarize, the lattice QCD calculations of strong-interaction processes that would 
be of fundamental importance for hypernuclear physics are: 

• An extraction of the scattering length and effective range for AN scattering from 
finite- volume simulations; 

• A determination of matrix elements of the axial current for E — > E, g-ss] 

• A determination of the r)NN and 77AA hairpin couplings in PQQCD; 
and of weak-interaction processes are: 

• A measurement of the coefficients of the weak nonleptonic two-point functions, and 

• A measurement of both the S— and P— wave amplitudes for the weak nonleptonic 
decays, e.g. A pii; 

• A measurement of the amplitude for the weak scattering AA^ NN. 

A lattice QCD determination of any one of these quantities would significantly improve 
our understanding of hyperons, and hypernuclei in general. 
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V. APPENDIX I : LUSCHER'S RELATIONS FROM EFT(7f) 



From the discussions in Refs. 2^, 13 it is clear that in the pionless theory describing non- 
relativistic baryons, each of mass M, the exact two-body elastic scattering amplitude in the 
continuum is 

where the C2n are the coefficients of operators with 2n derivatives acting on the nucleon 
fields (or equivalently with n time derivatives). Applying the PDS scheme to Iq gives 

ir'^ = {^^ + ^p) + 0{D-4) , (34) 



where p = V ME and hence 



^ ^ (35) 

M p cot — ip 



where the subtraction-scale dependence of the one-loop diagrams is exactly compensated by 
the corresponding dependence of the coefficients C2n(/i)- 5 is the energy-dependent S"— wave 
phase shift (we will only consider S*— wave scattering but this construction generalizes to all 
partial waves). 

We are interested in the energy-eigenvalues of this system placed in a box with sides of 
length L with periodic boundary conditions. We can find the energy-eigenvalues by requiring 
that the real part of the inverse scattering amplitude computed in the box vanishes, 

Re{ir'\L)) = , (36) 



where 

loiL) = ^^Ett^ ' (37) 

and the sum is over momenta, k, allowed on the lattice. 

The PDS value of this linearly-divergent integral is found by adding and subtracting 
the continuum limit of the integral evaluated at E = 0. One of the continuum integrals 
is evaluated with a momentum cut-off that is equal to the mode cut-off of the discrete 
summation, while the other is evaluated with dimensional- regularization and PDS. The 
subtraction-scheme dependence vanishes as expected, and we find that eq. (jH^ becomes 

j IJ I V 27r ' 

where the limit A^ — > oo is understood, and where A^ is the magnitude of the integer cut-off, 
related to the momentum cut-off in the box defined via |kmax. | = 2nh. jL. Solving for the 
values of p for which eq. (jHHj) is satisfied recovers Liischers result (s) 0, 2^. Given that 
Aj is set by the edge of the first Brillouin zone, clearly one must not be too cavalier about 
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taking the Aj oo limit in present day calculations however, in this work we will take 
the limit Aj oo. This derivation of Liischer's formula j27, makes explicit the analytic 
continuation of the zeta-functions that appear in his expressions. 
Using the effective-range expansion of p cot 6 



1 1 2 

a 2 



+ 



J9COt 5 

and assuming that L ^ a,r one finds that the lowest-energy eigenvalue is 

Ana 



En 



1 - Ci- + C2 



+ 



(39) 



(40) 



where the coefficients Ci 2 are 



Cl 



TT 



-2.837297 = -Zoo(l,0) 



TT 



C2 = c?-^Ehi = 6.375183 = ^ ( (Zoo(l, 0))' - ^oo(2, 0) ) 
^ i=^n IJI ^ 



(41) 



where the 2'oo(n, k) are defined in Refs. |27l. l2 

/ A. 

Zooin,k) = lirn^ ( ^ 

j,|j|7^fc 



- 5"^ 47rA,- 



(42) 



for integers n > 1 and fc. It is important to note that there are sign differences between 
our result in eq. (jlUj) and the analogous expression in Ref. ji^, H^- This is due to the 
opposite sign convention for the scattering length a, for which we use that common to 
nuclear physicists. 

We also need the energy of the next state in order to extract both the scattering length and 
the effective range. The construction is somewhat more complicated than that for the lowest 
state. One looks at the energy of the next-most energetic level that is invariant under the 
action of the complete cubic group; the space of such states is the Ai sector [27,,i28i]. It is easy 
to show from the decomposition of the full rotation group down to the cubic group [H^ that 
the energy shift of such states due to the two-body interaction receives contributions from 
/ = 0,4,6,.. interactions, where / is the angular momentum of the interaction. Therefore, 
to the order we are working in the momentum expansion, the shift in the energy-levels is 
sensitive to S*— wave interactions only. It is straightforward to show that the energy of the 
next-highest state in the Ai representation of the cubic group is 27. 28^ 



An^ 12tan5o / . , / , . 



.o + 4tan2 5o + ... + 0{L-'') 



This issue is discussed in Ref. 55]. We thank WiUiam Detmold and David Lin for discussions regarding 
this point. 

It is interesting to note that one can also use this expression to find the energy-eigenvalues in the limit 
where r <C i «C a, as may be appropriate for two- or few-body nuclear physics observables. We will 
explore this limit in forthcoming work. 
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1 



2tx 
1 



2 ^oo(i) 1) ; 



27r2 



TTCi + TT^C? - TT^Ca - 13 + \;U(\;\2 _ 



-0.061367 



^( (Zoo(l,l))2-6Zoo(2,l) 



47r4 



27r4 



vr^cs - TT^c? + 5 - ^ - 

j,|j|7^0,l 



2|j| 


2-1 


|j 




|2_ 1)2 



-0.354156 



(43) 



where 5q is the 5*— wave phase shift evaluated at the un per turbed lattice momentum, |p| = 
27r/L, and coefficients c'^ 2 ^^^^ those found by Liischer ji^, H^. We have given expressions 
for di 2 in terms of Ci^2 that can be straightforwardly evaluated and converge rapidly. 

To generalize the expressions in eq. ()40p and eq. ()43|) to AA^ scattering we make the 
replacement M and understand that 6q is the phase shift for AA^ scattering with 

the scattering lengths and effective ranges discussed in the text evaluated at the unperturbed 
momentum |p| = 2n/L. 



[1] For a review and comparative study of the different reaction mechanisms, see W. Alberico 

and G. Garbarino, Phys. Kept. 369, 1 (2002), nucl-th/01 12036, 
[2] J.J. Szymanski et al, Phys. Rev. C43, 849 (1991); H. Noumi et al, Phys. Rev. C52, 2936 

(1995); H. Outa et al., talk presented at the Vlllth Intl. Conf. on Hypernuclear and Strange 

Particle Physics (HYPOS), Jefferson Lab, Newport News, Virginia, Oct. 14-18, 2003, To be 

published in Nucl. Phys. A. 
[3] J. Golak, K. Miyagawa, H. Kamada, H. Witala, W. Glockle, A. Parreno, A. Ramos, and 

C. Bennhold, Phys. Rev. C55, 2196 (1997); ibid. Phys. Rev. C56, 2892 (1997). 
[4] A. Nogga, H. Kamada, and W. Glockle, Phys. Rev. Lett. 88, 172501 (2002), .iiucl-th/0112060| 
[5] E. Hiyama, Y. Kino, M. Kamimura, Prog. Part. Nucl. Phys. 51, 223 (2003). 
[6] C.L. Korpa, A.E.L. Dieperink and R.G.E. Timmermans, Phys. Rev. C65, 015208 (2001), 

nucl-th/0109072 

[7] J.F. Dubach, G.B. Feldman and B.R. Holstein, Annals Phys. 249, 146 (1996), 
tucl-th/9606003 

[8] A. Parreno, A. Ramos, and C. Bennhold, Phys. Rev. C56, 339 (1997) EucT-tIi/96 1 lOSOl 
A. Parreno and A. Ramos, Phys. Rev. C65 (2002) 015204 |imcl^Ii70104080[ 

[9] T. Inoue, K. Sasaki, M. Oka, Nucl. Phys. A670, 301 (2000). 
[10] A. Parreno, C. Bennhold and B.R. Holstein, nucl-th/0308074', 
[11] A. Parreno, C. Bennhold and B.R. Holstein, nucl-th/0308056,, 

[12] M. Fukugita, Y. Kuramashi, M. Okawa, H. Mino and A. Ukawa, Phys. Rev. Lett. 73, 2176 

(1994), |hep-lat/94076l2| Phys. Rev. D52, 3003 (1995), |hipTi[t79 501024, 
[13] http : / /nn-online . sci . kun . nl/index . html 

[14] B. Sechi-Zorn, B. Kehoe, J. Twitty and R.A. Burnstein, Phys. Rev. 175, 1735 (1968). 

[15] G. Alexander, U. Karshon, A. Shapira, G. Yekutieli, R. Engelmann, H. Filthuth and 

W. Lughofer, Phys. Rev. 173, 1452 (1968). 
[16] See, for instance, Y.C. Tang, in Proc. of the Int. Conf. on Hypernuclear Physics, Argonne Nat. 



22 



Lab. (Ed. A.R. Bodmer and L.G. Hyman), p. 276; R.C. Herndon and Y.C. Tang, Phys. Rev. 

159, 835 (1967); R.H. Dalitz, R.C. Herndon and Y.C. Tang, Nucl. Phys. B47, 109 (1972); 

T.H. Tan, Phys. Rev. Lett. 23, 395 (1969); J.J. de Swart et al, Springer Tracts in Modern 

Physics 60, 138 (1971). 
[17] S. Weinberg, Phys. Lett. B251, 288 (1990). 
[18] S. Weinberg, Nucl. Phys. B363, 3 (1991). 
[19] C. Ordonez and U. van Kolck, Phys. Lett. B291, 459 (1992). 
[20] C. Ordonez, L. Ray and U. van Kolck, Phys. Rev. Lett. 72, 1982 (1994). 
[21] D.B. Kaplan, M.J. Savage and M.B. Wise, Nucl. Phys. B478, 629 ( 1996) , |nuc]Tth/96 05002| 
[22] D.B. Kaplan, M.J. Savage and M.B. Wise, Phys. Lett. B424, 390 (1998), Eiuci-th/9801034l 
[23] D.B. Kaplan, M.J. Savage and M.B. Wise, Nucl. Phys. B534, 329 (1998) , Euci-th/9802075l 
[24] D.B. Kaplan, M.J. Savage and M.B. Wise, Phys. Rev. C59, 617 (1999), nucl-th/9804032, 
[25] S.R. Beane, RF. Bedaque, M.J. Savage and U. van Kolck, Nucl. Phys. A700, 377 (2002), 

nucl-th/0104030 

[26] M.C. Birse, J.A. McGovern and K.G. Richardson, Phys. Lett. B464, 169 (1999), 

|liep-ph/9807302 
[27] M. Liischer, Commun. Math. Phys. 105 153 (1986). 
[28] M. Liischer, Nucl. Phys. B354, 531 (1991). 

[29] H.W. Hammer, Nucl. Phys. A705, 173 (2002), nucl-th/011003ll 

[30] J.W. Chen and M.J. Savage, Phys. Rev. D65, 094001 (2002), |hep'-lat/0111050 

[31] S.R. Beane and M.J. Savage, Nucl. Phys. A709, 319 (2002), 1hep'-lat/0203003| 

[32] S.R. Beane and M.J. Savage, Nucl. Phys. B636, 291 (2002), hep-lat/0203028, 

[33] S.R. Beane and M.J. Savage, Phys. Lett. B535, 177 (2002), hep-lat/0202013, 

[34] S.R. Beane and M.J. Savage, Phys. Rev. D67, 054502 (2003), hep-l at/0210046 , 

[35] D. Arndt, S.R. Beane and M.J. Savage, Nucl. Phys. A726, 339 (2003), |nuc]Tth/ 0304004| 

[36] P.L. Cho, Nucl. Phys. B396, 183 (1993); [Erratum-ibid. B421, 683 (1994)], |hip^ 92Q8244| 

[37] R.F. Dashen, E. Jenkins and A.V. Manohar, Phys. Rev. D49, 4713 (1994), [Erratum-ibid. 

D51, 2489 (1995)], hep-ph/9310379, 
[38] D.B. Kaplan and M.J. Savage, Phys. Lett. B365, 244 (1996), |hep'-ph /950937Tl 
[39] S.R. Sharpe and N. Shoresh, Phys. Rev. D64, 114510 (2001), |hi^at/0108003l 
[40] L. Maiani and M. Testa, Phys. Lett. B245, 585 (1990). 
[41] K. Huang and C.N. Yang, Phys. Rev. 105, 767 (1957). 

[42] J.W. Chen, G. Rupak and M.J. Savage, Nucl. Phys. A653, 386 ( 1999) , Eucl-th/9902056l 
[43] U. van Kolck, Nucl. Phys. A645, 273 (1999), nucl-th/9808007, 

[44] P. van Baal, In the Boris loffe Festschrift, At the frontier of particle physics, vol. 2, 683-760, 

ed. by M. Shifman, World Scientific, | hep-ph /0008206| 
[45] J.F. Donoghue, E. Golowich and B.R. Holstein, Cambridge Monogr. Part. Phys. Nucl. Phys. 

Cosmol. 2, 1 (1992). 
[46] E. Jenkins, Nucl. Phys. B375, 561 (1992). 

[47] E. Jenkins and A.V. Manohar, UCSD-PTH-91-30, Talk presented at the Workshop on Effective 

Field Theories of the Standard Model, Dohogoko, Hungary, Aug 1991. 
[48] M.J. Savage and J. Walden, Phys. Rev. D55, 5376 (1997), |hep-ph/96112ld 
[49] R.P. Springer, "hep^ph/9508324 
[50] R.P. Springer, Phys. Lett. B461 (1999) 167. 

[51] A. Abd El-Hady and J. Tandean, Phys. Rev. D61, 114014 (2000), hep-ph/9908498l 

[52] J. Dai, R.F. Dashen, E. Jenkins and A.V. Manohar, Phys. Rev. D53, 273 (1996), 



23 



[hep-ph/9506273 

[53] R. Flores-Mendieta, E. Jenkins and A.V. Manohar, Phys. Rev. D58, 094028 (1998), 
hep-ph/9805416 

[54] L. Lellouch and M. Liischer, Commun. Math. Phys. 219, 31 (2001), hep-lat/0003023, 
[55] C.J. Lin, G. Martinelli, C.T. Sachrajda and M. Testa, Nucl. Phys. B619, 467 (2001), 
fhep-lat/0i64666| 

[56] J.E. Mandula, G. Zweig and J. Govaerts, Nucl. Phys. B228, 91 (1983). 



24 



